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ABSTRACT 

Generalizations of the Randall- Sundrum model containing a bulk scalar field 



O 

m 
o 

a: 

$ interacting with the curvature R through the general coupling are con- 

sidered. We derive the general form of the effective 4D potential for the spin-zero 
fields and show that in the mass matrix the radion mixes with the Kaluza-Klein 
^ ■ modes of the bulk scalar fluctuations. We demonstrate that it is possible to choose 

a non-trivial background form <&o(y) (where y is the extra dimension coordinate) for 
the bulk scalar field such that the exact Randall- Sundrum metric is preserved (i.e. 
such that there is no back-reaction). We compute the mass matrix for the radion 
and the KK modes of the excitations of the bulk scalar relative to the background 
configuration <&o(y) and find that the resulting mass matrix implies a non-zero 
value for the mass of the radion (identified as the state with the lowest eigenvalue 
of the scalar mass matrix). We find that this mass is suppressed relative to the 
Planck scale by the standard warp factor needed to explain the hierarchy puzzle, 
implying that a mass ~ 1 TeV is a natural order of magnitude for the radion mass. 
The general considerations are illustrated in the case of a model containing an i?$ 2 
interaction term. 
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1 Introduction 



Although the Standard Model (SM) of electroweak interactions describes success- 
fully almost all existing experimental data, the model suffers from many theoretical 
drawbacks. One of these is the hierarchy problem: namely, the SM cannot consis- 
tently accommodate the weak energy scale 0(1 TeV) and a much higher scale such 
as the Planck mass scale O(10 18 GeV). Therefore, it is widely accepted that the SM 
is only an effective low-energy theory embedded in some more fundamental high- 
scale theory that presumably could contain gravitational interactions. Recently 
many models that incorporate gravity have been proposed in the context of higher 
dimensional space-time. These models have received tremendous attention since 
they might provide a solution to the hierarchy problem. One of the most attractive 
attempts has been formulated by Randall and Sundrum [1] who postulated a 5D 
universe with two 4D surfaces ( "3-branes" ) with the following action: 

S = -Jd*x £ d y j ^ ( A + a) + ]T J\^\A k 5(y - y k )\ , (1) 

where R is the Ricci scalar, k 2 = 8ttG^ with G$ the Newton constant in 5D 
and A, Ai = Khid and A2 = A V i S are the cosmological constants in the bulk, on 
the hidden and visible branes, respectively. In the above, = 0,1,2,3,4) 

is the bulk metric and (gi)^ = (ghid)^u(x) = g^{x,y = yi = 0) and (g 2 )^ = 
(gvis)^u(x) = g^ u (x, y = y 2 = 1/2) (/i, v — 0,1, 2, 3) are the induced metrics on the 
branes. 

It turns out that if the bulk and brane cosmological constants are related by 



A = —, Khid = —^vis = — j- (2) 



and if periodic boundary conditions (y — > y + 1) with identification of (x,y) and 
(x, —y) are imposed, then the following metric is a solution of the 5D Einstein 
equations: 

ds 2 = e-^rj^dx^dx" - h\dy 2 , (3) 

where a(y) = m b [y(29(y) — 1) — 2(y — l/2)9(y — 1/2)]; b is a constant parame- 
ter that is not determined by the equations of motion. 
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Within the Randall- Sundrum (RS) model all the SM particles as well as the non- 
gravitational forces are assumed to be present on one of the 3-branes, the "visible 
brane" . Gravity lives on the visible brane, on the second brane (the "hidden brane" ) 
and in the bulk. All mass scales in the 5D theory are of the order of the Planck 
mass. By placing the SM fields on the visible brane, the initial 5D electroweak 
mass scale 0(M Pl ) is rescaled by an exponential suppression factor (the "warp 
factor") 7 = e - m obo/ 2 ; cl own to the weak scale 0(1 TeV) without any severe fine 
tuning. To achieve the necessary suppression, one needs m b /2 ~ 37. This is a 
great improvement compared to the original problem of accommodating both the 
weak and the Planck scale within a single theory. 

2 The radion effective potential 

A drawback of the RS model is the presence of a massless degree of freedom called 
the radion. There have been several attempts (see Refs. [2, 3, 4]) in the literature to 
generate the radion mass by introducing a bulk scalar field $ that would induce an 
appropriate radion potential. Here we will derive the general form of the potential 
within a class of models containing the bulk scalar interacting with gravity in the 
following manner: 



where we have introduced the bulk potential V(<&) and the brane potentials Vi($) = 
Vhid(^) and V 2 ($) = V vis (&). In addition to the standard scalar kinetic-energy 
term, we have allowed for a general coupling of the bulk scalar to gravity through 
the aRf(&) interaction term. Since we would like to preserve the explanation of 
the hierarchy proposed by Randall and Sundrum, we will also require that the RS 
metric (3) remain an exact solution 1 of the Einstein equations even in the presence 
of the bulk scalar $. Therefore, it is useful to separate out in the action (4) both 

1 One can, of course, consider slight modifications of the RS metric that would also solve the 
hierarchy problem in a similar manner. However, in this paper we would like to discuss a scenario 
with exactly the same metric as in the original RS model. 



s = 



J d A x j^dy^[\g 



-_ R 





(4) 
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the bulk (A) and brane (A^d, A^ s ) cosmological constants that satisfy the RS 
conditions, Eq. (2). 

In order to identify the radion, it is sufficient to consider scalar excitations of the 
metric around the background RS solutions. Hereafter, we will adopt the following 
parameterization (see Refs. [5, 6]) of the metric fluctuations: 



ds 2 = e -2* {y) -2 b{x) e^ (fW + y))dx , dx u -bl(i + 2b(x)e 2 ^) 2 dy 2 , (5) 



where h^x^y) and b(x) are related to the graviton 2 and the radion, respectively. 
Then from — ^J\g\R/ (2k 2 ) in the action (4) one obtains [after expanding in powers 
of b(x)} the kinetic term for the radion: 3 



It is easy to verify that if interactions of the scalar field <3> are switched off, then 
there is no potential for b(x) and consequently the radion would be massless. 

The bulk scalar has been introduced here in order to generate a non-trivial 
potential for the radion. However, in general the presence of the scalar leads to a 
non-trivial interaction potential between the radion and the scalar in addition to 
the appearance of a radion potential. Therefore, the strategy that we will follow 
here will be to determine the background scalar configuration $(?/) such that the 
RS background metric is preserved and then to expand the action (4) around it. 
First, one has to solve the Einstein equations together with the equation of motion 
for $. Let us start with the Einstein equations, keeping in mind that we would like 
to preserve the RS metric as a vacuum solution. We write 



where Gij is the Einstein tensor, denotes the RS contribution to the energy- 

momentum tensor and [ST)ij contains all new contributions emerging from inter- 
actions of the scalar $. It is useful to first calculate all the extra (compared to the 

2 In the following we will not discuss interactions with gravitons as they will not influence the 
potential for scalar degrees of freedom. 

3 Hereafter, the flat metric r\ liv will be assumed whenever repeated indices are summed. 





Gij — K 



2 \ T (RS) 



(7) 



4 



pure RS model) contributions to the energy-momentum tensor. It is easy to show 
that: 



(«r)«=TW + 2a{Dy[/($)]-Gy/($)} + f £ V k m 9k ), u 5^8(y-y k ), (8) 



1 



where 



Tip = V^V^-ft, 



& o fc= i : 



= ViVjX-gijcfVkVtX. 
For the RS background metric we obtain: 

-2„-2cr/o '2 



(9) 
(10) 



G 



36 - 2 e- 2CT (2a 2 - ( 7")^ 
6a' 2 



and 



Damn 



b 2 e- 2CT (-3a'/($)' + /(*)>, 







(12) 



4a ft*) ) ' 

where here, and in what follows, the prime denotes differentiation with respect to 
the 5th coordinate, y. 

Since we demand that the RS metric be preserved even when the scalar is 
present (no back- reaction from the scalar), we have to require that the extra con- 
tributions to the energy momentum tensor (calculated using the RS metric) vanish: 



(«%($) =0. 



(13) 



Since we want to find a background solution for $ that satisfies 4D Lorentz invari- 
ance, we will assume that the solution is only a function of the extra dimension 
coordinate, y. The (//, v) and (5, 5) components of Eq. (13) read, respectively: 



:$') 2 + 12a(2o-' 2 - a")f - Ylaa f + 4af"+ 



2bl 



V($) + - £ V k (*)6(y-y k ) 



J k =h2 



$') 2 - 24«(j' 2 / + lQaaf - 26 2 K($) 








(14) 
(15) 



Note that since $(y) should be a continuous function, the above equations imply 
$'(y) = and V(<&) = V vis (&) = V hid ($) = for a = 0. Therefore, introduction 
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of the extra coupling aRf(&) is essential in order to obtain a no-back reaction 
solution, (8T)ij($) = 0. 

In addition, $ must satisfy the following equation of motion: 



+ 4a' & + 4a(5cx 2 - 2a")^- + b 2 



0, (16) 



where the RS metric was used. 4 Once the vacuum solution ($ ) is determined, we 
expand the action, Eq. (4), adopting the parameterization of the scalar fluctuations 
of the metric given in Eq. (5) and the following definition for the <3> quantum 
fluctuation: 

$0,y) ^%(y) + <f)(x,y). (17) 

Then, in order to determine the effective 4D potential for the scalar degrees of 
freedom, we collect all non-derivative contributions to the d 4 x integrand in the 
action of Eq. (4) containing b(x) and <p(x,y). 

In other words, we expand the theory defined by the action (4) around the 
vacuum solutions for the metric (the RS solution) and the bulk scalar field $ [the 
solution of Eqs. (13) and (16)] in terms of the scalar fluctuation of the metric, b(x), 
and fluctuation of the scalar field, (f)(x,y). 

First, let us calculate the Ricci scalar for the metric (5) and collect all the terms 
containing derivatives with respect the the extra component y: 

D 20 , 2 8 a" 

R -¥ a ~¥ TT2b^ + "' (18) 

where ellipses contain only (x, y)-derivatives of the graviton and re-derivatives of 
the radion. Since we are going to calculate the potential, derivatives of b(x) will be 
dropped hereafter. As has already been mentioned, we will not consider fluctuations 
of the rjnv part of the metric. Therefore, we will also neglect all terms containing 
h^x^y). 5 

4 Both the Einstein equations (13) and the equation of motion (16) constrain the scalar $. 
However, it can be verified that the equations are not independent; a certain linear combination 
of derivatives with respect to y of the (/i, v) and (5, 5) components of Eq. (13) is proportional to 
Eq. (16). Since we have not specified the potential V($>), we can find a consistent solution both 
for $ and the potential. Details of the derivations for /($) = 3/32 <I> 2 will be presented in Sec. 3. 

5 Adopting the traceless gauge, h£ = 0, one can eliminate possible mixing between b or 4> 
and hp V . Consequently, graviton interactions cannot influence scalar masses. Thus, h^ v will be 
suppressed in the following. 
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Using the contributions to the Ricci scalar displayed in Eq. (18), one gets the 
following form of the effective 4D potential from the action (4): 



dy e - 4 - 4be2CT x 



2 |_&o(l + 2be 2 °) 



:$ o + 0) 2 + 26 o (l + 26e iCT )\/(<l>o + 0) 



+ 



K 2 b 



5(1 + 2be 2a )a' 2 -2a" 1 + 2/t 2 a/($ + 0) 



+ 



A6o(l + 26e 2CT ) + £ [A fc + V k (Q + <P)}S(y - y k ) 



(19) 



fc=l,2 



where <3> = $o(y) denotes the vacuum solution (that preserves the RS metric) for 
the scalar $. Note that $ is determined as a solution of the equations of motion 
for the RS background metric. As a result, it does not contain any dependence 
on the radion field b(x). It is easy to verify that contributions from the pure RS 
model to V e s(b, <fi) vanish when the relations (2) are satisfied: a non-trivial potential 
requires an extension of the minimal RS model. 

Next, it is easy to show from Eq. (8) that if $ is independent of x then the 
following identity holds: 



-($') 2 + 26 \/($) 



24a / / o ;/\ „ . T . 
+ (2a 2 - a ) /(<£>) 



bo ibo 

(-to'W)' + /(<&)") + ^ E ^(*)% - Vk) • (20) 

°0 °0 fc=1>2 

Multiplying the above equation by exp{— 4cx}(&o/4), integrating the (— 3er'/ + /") 
term by parts and using the RS relations, Eq. (2), one obtains the following simple 
relation between (5T)^($ ) and the minimum of the effective potential of Eq. (19) 
at its [b(x) = 0, <f>(x, y) = 0] minimum: 



M0,0) = 7 r\dye-^{5TT^ ). 



(21) 



Note that in Eq. (21) we employ the trace from Eq. (20) calculated for the back- 
ground solution $o- Since the no-back- reaction requirement, Eq. (13), implies 
(5T)^($ ) = 0, the relation (21) shows that the effective potential must vanish at 
the minimum 



Kff(0,0) =0. 



(22) 
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It is straightforward to verify that linear terms in b(x) and <f>(x, y) disappear by 
virtue of Eq. (13) and Eq. (16), respectively. 

In order to determine scalar masses one has to expand the action (4) up to terms 
quadratic in b and 0. First, let us define the K-K modes of the scalar fluctuations: 

J n (y) 



X) 



b\ 12 



(23) 



with orthonormal functions J n {y): 

dye- 2 ^Uy)J m (y)=5 n 



(24) 



The resulting mass terms are the following: 



\ ( r Vm) 



Ml {M^) n \ ( r 
(M^) m (M 2 ) mn J \ <p r . 



(25) 



where r is the canonically normalized radion [see Eq. (6)]: 

6 X 1 / 2 



r(x) 



K 2 m 



D m b 



' 2 b{x) 



(26) 



Inputing the equation of motion (16), the elements of the mass matrix read: 



Mi = 



^2^0 e -m 6o 



3 b 



\ (*o) 2 dy 



(27) 



$o + 2a + 20a O -' 2 ^($ ) + &o^($o) 



(28) 



(M 2 ) mn 



h 2 I 

°0 J - 2 



[e-^\j n ( y )J m (y)x 



4a (5a 2 - 2a" 



J ' n {y) J 'm(y) 



d 2 f 
) — - 

d$ 2 



d 2 V 1 



d 2 V k 



d§ 2 b fc tr 2 d ® 2 



5{y-yk) 



+ 



(29) 



Before we can estimate the size of the elements of the mass matrix, we must 
discuss first the constraint that is imposed on the model by the requirement of 
maintaining the standard strength of classical 4D gravity. Adopting the metric 
defined by Eq. (5), one can calculate the coefficient of the 4D Ricci scalar obtained 
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for g^y = rj^+h^. In order to reproduce the standard result, the coefficient should 
be M Pl /2. The resulting constraint is: 



where M P \ ~ 2 x 10 18 GeV is the reduced Planck mass and 7 = exp(— m &o/2). 
In order to solve the hierarchy problem, one needs mobo/2 ~ 37. Therefore, terms 
of order 7 2 can be safely neglected in Eq. (30). It is clear that the most natural 
scenario 6 emerges when all the mass parameters of the 5D theory are of the order 
of Mpi. In this case, the elements of the scalar mass matrix defined by Eq. (25) 
are of the following order of magnitude: 



where a and b are calculable coefficients of the order of 1. It is clear that for 
m b /2 ~ 37 the lowest scalar mass is of order: 



There are two essential conclusions. First, we see that the lowest scalar mass 
receives the standard suppression from the warp factor 7 = exp(— mo&o/2) that is 
necessary for the solution of the hierarchy problem. As a result, the mass expected 
in the presence of a bulk scalar Rf(&) interaction is of the order of the electroweak 
scale. Second, it is clear that in order to find precise values for scalar masses 
for any particular choice of the interaction function /($), one has to take into 
account the mixing between the radion and the Kaluza-Klein modes of the bulk 
scalar fluctuation. To fully explore the phenomenology of the theory, it would be 
necessary to calculate all the entries of the mass matrix; however, this is beyond 
the scope of this paper. 

Finally, we close this section by reiterating the fact that if there is no Rf(&) 
interaction, i.e. if a — 0, then necessarily <&' (y) = and V(Q ) = V vis (& ) = 

Vhidi&o) — 0, which in turn would lead to a vanishing mass matrix. 

6 Of course, an appropriate cancellation between contributions coming from parameters that 
differ even by many orders of magnitude is in principle also possible. However, since we would 
like to preserve the solution of the hierarchy problem proposed by Randall and Sundrum, we 
should assume that all the mass parameters in the fundamental 5D theory are of the same order. 
Then, the only necessary fine tuning is to keep mobo/2 <~ 37. 




(30) 




(31) 



-fMpi(a-b) 1 / 2 246 GeV 
v'2 " ^ v '2 



(32) 
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3 The R<$> 2 interaction 



In this section, we will illustrate the general discussion from Section 2, choosing a 
specific form of the interaction between the bulk scalar and gravity: 

3 



/(*) = i* 2 



(33) 



The function /($) has been normalized such that a — — 1 corresponds to a 5D 
conformally invariant interaction. This coupling was discussed in various contexts 
by many authors in the past, see e.g. Ref. [7]. 

In this case, the conditions for (5T)ij(&) = 0, Eq. (13), as written out in 
Eqs. (14) and (15), read: 



-3a'($ 2 )' + ($ 2 )" + ^($') 2 + 3(2a' 2 - cx> 2 



166 2 
3a 



1 



V(*) + r £ V k (*)6(y-y k ) 

°0 k =l,2 



0; 
0. 



(34) 
(35) 



Eliminating V(<&) from Eqs. (34,35), one obtains the following equation for $: 

«/($*)' + ($ 2 )" + i^($') 2 - 3a" V + ^ E -Vk) = 0- (36) 



Away from the branes, we find the solution: 

$o(y) =d\l- ce-^ 



(37) 



where c, d are integration constants and (3 = 2 + 8/ (3a) is required for consistency. 
If is not an integer, we must also demand that 1 — ce~ a ^ > in order that 
$o(y) be well-defined. Recalling that a" = 2m b [5(y) — S(y — 1/2)], and noting 
that ($ 2 )" will contain a term proportional to a", the conditions that the coefficients 
of the (^-functions in Eq. (36) vanish reduce to: 



m o d 2 g(c)+P(p-2)V hid (0) = 
m d 2 g( C1 )-P((3-2)V vis (^) = 0, 



(38) 
(39) 
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where 7 = exp(-m 6 /2), g(x) = (l-x)^ 1 [x(2 + 3(3) - 3/3], V hid (0) = V hid [<S> (0)} 
and V vis (j^J = V V i S $0 (|) , and we have introduced the notation V V i S = V2 and 
V hid = V 1 . 

Insertion of the solution <& (y) into, for example, Eq. (35) fixes the form of the 
bulk potential: 



V(*) 



3 2^2 J 4 
-mi§ z < — 



~d 



+ 



0-2 

6(3 2 



~d 



(40) 



In addition, <3> must satisfy its equation of motion as obtained from Eq. (16) for 
the form Eq. (33): 



4a $ - $ + — (5a 2 - 2a )$ + 6 2 



0. (41) 



It is easily verified that the bulk form for &o(y), Eq. (36), also satisfies the equation 
of motion in the bulk, Eq. (41). However, cancellation of the 5(y) and S(y — 
1/2) brane delta function pieces yields matching conditions that are different from 
Eqs. (38,39). For consistency of Eq. (41) we require: 



2m dh{c)+(3{(3-2)V M M = 



0, 



2m dh( C1 )-f3(f3-2)V vis [- 



(42) 
(43) 



where h(x) = (1 -x)^ 1 [ a ;(2 + 3/7) -4/3], V Md (0) = dV hld /d^ = „ m and V' ms (±) = 
dV ms /d^ = ^_y 

Equations (38) and (42) can be solved for the parameter c in terms of Vhi d (0) 
and Vfr id (0). Two solutions are possible for c, specified by 



<k = fi(P,-Rhid) « = 1, 2 , 



where 



Rhid = 



(44) 



(45) 



4m V hid (0) ' 

The functions fi denote the two possible solutions of the quadratic equations for c: 

-B ± VB 2 - 4AC 



MP,R) 



2A 



(46) 
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where i — 1,2 corresponds to the + and — signs in front of the square root, 
respectively. The quantities A, B, C in the above are given by: 

2 + 3/5 



A = 
B = 
C = 



-(2 + 3/5) 
(2 + 3/3) 



1 + 



W-2) 



1 + 



R(P-2) 



3/3 



-P 



3 + 



16 



R{(3-2)_ 

Positivity of B 2 - AAC leads to the requirement 1 - (2 + 3(3) /R > 0. 

Once c is determined, we can compute 07 in terms of V' vis and V vis [^j using 
Eqs. (39) and (43). One finds 



Cili = fj (P, Rvis) for ? = 1,2 



where 



Rvis — 



V' 2 (±) 
v vis y 2 ) 



(47) 



(48) 



^m V vis (i) 

and the appropriate branch j = 1 or j = 2 is determined by the need to obtain 
a very small value for the warp factor 7j, i.e. fj{f3,R V is) ~ 0. The latter is 
most straightforwardly achieved by requiring C V i S ~ and choosing the j = 1, + 
(j = 2, — ) solution for fj(/3, R V i S ) [see Eq. (46)] for B vis > (B vis < 0), respectively. 
From Eq. (47), the requirement that C vis ~ is equivalent to 



or Rvis — 



(49) 



3(1 - \(3) ■ 

For i?m S as above, 7 the positivity of \J B 2 is — AA vis C V i S is automatic so long as B vis 
is not extremely tiny. For R vis as given in Eq. (49), one finds B vis ~ 3/3/2 — 1, so 
that we must use the j = 1, + solution for /3 > 2/3 and the j = 2, — solution for 
(3 < 2/3. For (3 = 2/3, the choice becomes ambiguous. 

For convenience, in what follows we denote by A vis , B vis , C V i S the values of 
A, B, C for R = R vis and by A hid , B hid , C hid the values of A, B, C for R = -R h id- 

7 For R V i S = , \ = r(f3) exactly, 7 = 0. As R V i S changes (for a fixed Rhid) from a value 
slightly larger than this to a value slightly smaller, 7 will change sign. Whether a positive or 
negative shift relative to r{[3) is required to give -fi > is determined by the sign of a, which can 
vary. 
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From Eqs. (44) and (47), the warp factor 7$ (and hence the distance between 
the branes) for a given solution Cj is given by 

7 ,^ e -^= /'jf'**' for -1,2, (50) 

with j = 1 for /3 > 2/3 and j = 2 for (3 < 2/3, as discussed earlier. In practice, 
we will require that 7^ = 7 = e~ 37 for either choice of i Further, one can use (for 
example) Eq. (42) to determine df. 

P(J3-2)VU0) rM x 
at " 2moM/,(/3,- J R^)]- l ° J 

Once dj, q and 7$ have been fixed as specified above, Eqs. (38) and (39) imply a 
consistency constraint on the model parameters: 

V hid (0) _ g[fi((3,-R hid )} g[fi(P,-R hid )} 



V ms (f) = ~ 9[fj(P,Rvis)] ~ 3/5 



(52) 



where 7$ ~ has been used to obtain the last approximate form. Using 7, ~ 0, 
Eqs. (39) and (43) also simplify to 

K,fi)~-3?4, (53) 



,2/ /3-2' OTS V2/ /3-2' 

with dj as given in Eq. (51) for solution branch i. Note that, for a given value of 
(3 and choice of solution branch i, fixing V V i S (j^J corresponds to fixing the normal- 
ization di of $0 in terms of m and that fixing also V^ s (|) then fixes both m and 
di. 

Finally, it is important to note that the definition of Rhid, Eq. (45), yields the 
following constraint on the relative signs of R^d and Vhi d (0): 

V'^M = 4m R h i d V hid (0) > . (54) 

Using Eqs. (53) and (52), the condition Eq. (54) can be converted to a requirement 
expressed entirely in terms of Rhid and (3 for a given solution branch i: 

Rhidg(fi(P, —Rhid) n /__>. 
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The conformal limit of a — — 1 (/3 — —2/3) requires special treatment, 8 since 
for this choice A hid = 0. In this case, g(x) = 2(1 — x)~ 4 and h(x) = |(1 — x)~ 5 / 2 . 
Eqs. (38) and (42) then yield 

m d 2 (l - c)- 4 = ~V hid (0) , m d(l - c)- 5 / 2 = —V^M , (56) 

respectively, from which we conclude that V^ d (0) < and V^ d (0) < are required, 
which also implies that [see Eqs. (45) and (54)] R hid < 0. By combining Eq. (56) 
and Eq. (45) we find 

c=l + J-. (57) 
tihid 

Note that Rhid < —2 is required for < c < 1, but that c is negative for —2 < 
Rhid < 0. There is nothing obvious to forbid this latter choice since (1 - ce-^)) 
will automatically be positive for all y if c < 0. In an analogous spirit, utilizing 
Eqs. (39) and (43), one can show that 

try = 1 - -i- . (58) 



1 rv.s 



Combining Eqs. (57) and (58), one obtains the following result for the warp factor: 

Rhid Rvis 2 /rn\ 

7 = R .9 • ( 59 ) 

-fkiis J^hid ~r ^ 

In order to have a phenomenologically acceptable small value for the warp factor 
7, either R vis ~ 2 or it^ ~ is required. The remaining constraint [the analog of 
Eq.(52)] in this case reads: 



V hid (0) (R 



l-hid 



Rv 



4 



(60) 



vms y 2 J 

Combining this equation with the earlier-noted constraint that V^(0) < results 
in the requirement that V V i S (j^ > 0, implying that R vis > 0. Combined with 
Eq. (59) and the requirement that 7 > 0, the only allowed choices are: 

(0 < Rvis < 2 and - 2 < R hid < 0) or (R vis > 2 and R hid < -2) . (61) 



8 In addition to the general solutions discussed in the main text for this case, there exists a 
special background solution, $o(y) oc e 3 / 2<T ^\ for which there are no matching conditions since 
the solution satisfies all the necessary equations everywhere, including the boundaries. For this 
particular special conformally symmetric case, substituting the form of <f>0i as given above, into 
Eq. (35) leads to a vanishing bulk potential, V(&o) = 0. Similarly, Eq. (36) gives Vi.2( < f ) ) = 0. 
Because all the potentials are zero, one finds = 0. We are only interested in cases for which 
a non-zero mass is generated for the radion. 
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(Note that Eq. (49) does not apply for the conformal choice of a.) For R v 
as generally needed for small 7, one finds 

1\ 9 , „ _ 9 



V„. 



9 

-m d 2 , V hid (0) = ni 



m d R hid . 



2, 



(62) 



In Fig. 3, we display & (y)/d as a function of y for three cases. In all cases, we 
have chosen input parameters so that 9 m b /2 ~ 37 as required for the warp factor 
1 TeV/Mpi. In the first case, we have taken (3 ~ +2/3, equivalent 



7 = e" 



-m b /2 



to R vis = 4 from Eq. (49), and R hid = +1. For this choice, c ~ 0.7835. In the 
second (third) cases, we make the conformal choice of (3 — — 2/3, take R vis ~ 2 
(for small 7) and employ R hid = —4 (Rhid = —1) for which c = 1/2 (c = —1). In 
the (3 = 2/3 case (which is representative of cases with (3 > 0), we see that $o(y) 
is repulsed from the hidden brane located at y — 0. The second (third) case is 
representative of a (3 = —2/3 case for which $o(y) is strongly peaked on (strongly 
repulsed from) the hidden brane. 

A useful cross-check is to adopt the explicit form of the solution (37) to verify 
that indeed V(0,0) vanishes as predicted by Eq. (22). In order to calculate the 
radion potential at the minimum we use Eq. (35) to eliminate V(<&) in the general 
formula (19). The result is the simplified expression 



Veff (0,0) 



10mQ&o 

8m 
7^2 



dye^(%(y)Y 



$o(0) ; 



$0 



1 

b~o 

-2m 6o 



dye-^(^(y)Y+ (63) 



+ V hid (0) + V vis [ - 1 ( 



-2mo&o 



Then, inserting the solution Eq. (37) into Eq. (63), one obtains: 

(3((3-2)V eff (0,0) 

= 'mo d 2 (l - c)^- 1 [c(2 + 3/3) - 3/3] + /?(/? - 2)\4, d (0) _ 



-7 



m rf 2 (l - c^^" 1 [c 7 (2 + 3P) - 3(3} - (3((3 - 2)V ms (l/2) 



m d 2 g(c)+(3((3-2)V hid (0) 



-7 



m d 2 g( C1 ) - (3((3 -2)V vis (l/2) 



(64) 



9 We reemphasize that mo&o/2 is calculable in terms of the input parameters using Eq. (50) or 
Eq. (59). 
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Figure 1: The normalized scalar background function Q (y)/d 
mobo/2 ~ 37 in three cases: (top) Rua = 1 and (5 ~ +2/3 (equiva^ 
(middle) /3 = —2/3 and = —4 (implying c = +0.5); (bottom) (3 = —2/3 and 
-R/iid = —1 (implying c = —1). In the latter two (conformal) cases, m bo/2 = 37 
requires R V i S ~ 2. 
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The two bracketed expressions on the right hand side of the above equation vanish 
by virtue of the fact that they correspond precisely to the expressions appear- 
ing in the matching conditions that emerge from the requirement of (ST)^ U = 0, 
Eqs. (38,39). 

Finally, if we insert the solution Eq. (37) into the formula for the radion mass, 
Eq. (27), one obtains [after employing Eqs. (38) and (39)] the following result for 



Ml 



9 

■h/rl ^ 2 2 

M r = -k m 7 



-V hid (0) 



c{2-(3)+(3 



V m 



1\ c 7 (2-/3)+/3 " 
2 J c 7 (2 + 2,(5) - 3/3 



(65) 



c(2 + 3/3) - 3(3 

where c and cy are given by Eqs. (44) and (47) or, in the conformal case, Eqs. (57) 
and (58). Note the presence of the warp factor 7 that reduces the radion mass from 
the typical 5D scale [a natural choice for which would be O (Mpi) if V^O) and 
V V is (l) are O (M pl )] down to the electroweak scale. Of course, Eq. (27) guarantees 
that M 2 > for any given form of /($) and any background solution $o(?/) that 
is real and fully consistent. The conditions for reality and consistency are: (a) 
the constraint Eq. (52) or, in the conformal case, Eq. (60) is satisfied; (b) the 
positivity condition Eq. (54) is satisfied; and (c) 1 — c > when 1//5 is not an 
integer. However, substantial variation of M 2 is possible. In particular, M 2 = 
at special points when considered as a function of at fixed (3. A final form for 
M 2 can be derived in the 7 1 limit by using Eqs. (52) and (53) in Eq. (65): 

g(c) c(2-P)+P~ 



2 (\ 

M 2 r = -K 2 m 0l 2 V ms 

= 3(^7) 7J- 



1 - 



P c(2 + 3[3)-3[3 



:i-c 



4-1 



!-c + c- -1 



(66) 



This shows the fundamental importance of the scale of V v , 



in determining 



M r 2 . Of course, one should continue to keep in mind the relation Eq. (53) between 
V V i S {^j and m d 2 as well as the relation Eq. (52) between V v i s (|) and 14^(0). 

To illustrate how M 2 depends upon the parameters and how other constraints 
enter, consider first the conformal case with (3 = —2/3. For R vis ~ 2 (so that 
7< 1) and using Eqs. (57) and (60), one can rewrite Eq. (66) in terms of V vi 
(which, as noted earlier, fixes the value of d in terms of m ). The result is: 



M'l = ^m 0l 2 V vis ( ^ ) K(R hid ) = h^Km d) 2 K(R hld ) , 



(67) 
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where 

K(R hid ) = ^R{ ld + l -R\ id + 1 > . (68) 

In Fig. 2, we plot K(R hid ) as a function of Rhu- We observe that M 2 > ev- 
erywhere except at Rhi d = —2 (a point where c = is required by Eq. (57) and 
$o(y) becomes trivial). Our two earlier (3 = —2/3 plots of the wave function thus 
correspond to choices for which M 2 > 0. 




Figure 2: The function K(Rhi d ) is plotted as a function of Rma in the (allowed) 
R hi d < region. 

For fixed (3 ^ —2/3, and 7 <C 1, M 2 is approximately a function of R hid only, 
where Rud is to be restricted to those values such that a given solution c\ or c 2 sat- 
isfies the other consistency constraints. We explore the behavior of M 2 as follows. 
First, recall that Kis (|), ^^(0), R V i S , Rhid and (3 are the input model parameters. 
At fixed (3, in order to obtain 7 1, as desired on phenomenological grounds, we 
adjust R vis according to Eq. (49). Then, for a chosen R hid we calculate q using 
Eq.(44). In fact, there are the two solution branches, c 1:2 (Rhid) — fipiP, — Rhid)- 
Positivity of B\ id - AA hid C hid , as determined by the sign of 1 - (2 + 3f3)/(-R hi d), 
requires R hid < Min [-(2 + 3/3), 0] or R hid > Max[0, -(2 + 3/3)] (the limit cross 
over taking place at the conformal point of f3 — —2/3). Once C\ [c<i) is chosen, if 
1/(3 is not an integer we check to see if 1 — ci > (1 — c 2 > 0), as required for 
a well-defined $0(2/) = d(l — ce~ 2<j( ^) 1 / /3 . Finally, for any values of (3 and Rud, 
and given choice of branch i, we check the positivity requirement of Eq. (54). To 
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compute M 2 , we adopt Eq. (66), which takes into account the consistency con- 
straint (52) such that 14^(0) is expressed in terms of V vis (j^. Eq. (66) makes it 
clear that, for a given value of (3, M 2 is proportional to V vis (j^J [which fixes m df 
through Eq. (53)] and depends non-trivially on Rhid, the other parameters being 
held fixed. Remarkably, one finds (M 2 )j > (as expected) so long as: (a) Rhid is 
such that 1 + (2 + 3f3)/R hid > (so that Ci j2 are real); (b) 1 — q > when 1/(3 is 
not an integer; and (c) the positivity condition, Eq. (54), which we abbreviate as 
RhidVhid > 0, is satisfied. There are many different cases. Here we simply describe 
a couple of illustrative possibilities. 

K' 

400 
200 



-200 
-400 

Figure 3: For (3 = 2/3, we plot [K'] 2 = [M 2 /{§(/tm d7) 2 }] 2 , see Eq. (66), as a 
function of R hid in the R hid < -(2 + 3(3) = -4 region. Both 1 - c 2 and (RhidVhid) 2 
display the same sign changes and singular behavior as [M r 2 ] 2 . Only the [K'} 2 > 
region corresponds to a solution consistent with all constraints. 

Consider first two choices of (3 such that 2 + 3(3 > 0. 

• For (3 = 2/3, any value of R hid > gives 1 - c lj2 > and (RhidVud)i,2 > 
and substantial values for the corresponding (M 2 ) x and (M 2 ) 2 . For R hid < 
— (2 + 3(3) = —4, any value of Rhid gives 1 — c\ > and substantial (M r 2 )i. 
In all these cases, M 2 is a relatively smooth function of Rhid- 

However, in the Rhid < —4 region, 1— c 2 and (RhidVhid)2 are both only positive 
for —4.5 < Rhid < —4 and (M 2 ) 2 varies rapidly, as illustrated in Fig. 3. This 
case illustrates the extreme sensitivity that M 2 can have to Rhu and shows 
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that very large and very small values of M 2 are quite possible. 

• For (3 = 1/2, since 1/(3 is an integer, the sign of 1 — c is unconstrained. 
One finds (RhidVhid)i,2 > for all R^ d > and (M 2 ) lj2 > and behaves 
smoothly. For all Rud < —7/2, {RhidYhid)\ < and this solution branch 
fails the positivity constraint. In the small region —4.65 < Rhid < —7/2, 
{RhidVudh > and (M r 2 ) 2 > 0, blowing up at R hid 4.65. 

We also briefly describe one interesting 2 + 3/5 < case. 

• For (3 — — 1, again the 1 — c > constraint is not necessary. For all Rhid > 
-(2 + 3/?) = 1, both (RhidVhid)i < and {RhidVhidh < so that the c x and c 2 
solutions both fail the positivity constraint. In contrast, for all Rud < one 
finds {R h idVhid)i > and, of course, (M r 2 )i > 0. However, (R h i d V hid ) 2 > 
only for —0.32 < R hid < and R^d ^, —1.78 — in these two ranges (M 2 ) 2 > 
aside from a zero at the top end, Rud — —1-78, of the 2nd range. 

Overall it is clear that there is a large range of possible models that satisfy all the 
constraints necessary for exact Randall-Sundrum metric with positive radion mass- 
squared. (Some particular choices are somewhat more fine-tuned than others.) 
We have not understood how to choose between the various models; possibly, the 
conformal models with (3 = —2/3 should be regarded as the more attractive. 

4 Conclusions 

We have considered a class of generalizations of the Randall-Sundrum model con- 
taining a bulk scalar field $. We demonstrated that no-back reaction from the 
scalar on the Randall-Sundrum metric solution requires the existence of an ex- 
tra interaction between gravity and the scalar. Here, we considered the coupling 
i?/($). A general form of the potential for the fluctuation of the compactification 
volume (the radion) and the Kaluza-Klein modes of the excitation of the bulk scalar 
was derived and the mass matrix was determined. In order to obtain the values 
of scalar masses, one has to take into account the mixing between the radion and 
the Kaluza-Klein modes of the fluctuation of the bulk scalar. We demonstrated 
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that a non-zero mass for the lowest eigenstate (which we identify with the phys- 
ical radion) can be generated using a choice of the background bulk field, $o(y)> 
that preserves the RS metric (no back-reaction). We found that the radion mass 
receives the same suppression from the warp factor that is necessary to explain the 
hierarchy puzzle. Thus, ~ 1 TeV is a natural order of magnitude for the radion 
mass. Finally, we illustrated the general scenario for the case of /($) oc $ 2 , for 
which the scalar background solution that preserves the Randall-Sundrum metric 
was explicitly found. 

Since the mass-squared matrix for the radion and KK bulk scalar excitations 
is non-diagonal, it is clear that the introduction of Higgs-radion mixing on the 
visible brane through a term in the Lagrangian of the form ^ \g V i S \R(g vis )H^ H , 
as considered for example in [8], would result in a complicated situation where the 
Higgs field, the radion and the KK excitations of the bulk scalar would all mix. 
A phenomenological study of the magnitudes of the various mixings, as a function 
of the available parameters, would be required to understand the extent to which 
the phenomenology of the physical Higgs boson eigenstate can be modified. Such 
a study is beyond the scope of this paper. 
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